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Abstract 

A microscopic model of indirect exchange interaction between transition metal impurities 
in dilute magnetic semiconductors (DMS) is proposed. The hybridization of the impurity d- 
electrons with the heavy hole band states is largely responsible for the transfer of electrons 
between the impurities, whereas Hund rule for the electron occupation of the impurity d-shells 
makes the transfer spin selective. The model is applied to such systems as n— type GaN:Mn and 
p— type (Ga,Mn)As, p— type (Ga,Mn)P. In n— type DMS with Mn^"^/^'*' impurities the exchange 
mechanisms is rather close to the kinematic exchange proposed by Zener for mixed-valence Mn 
ions. In p— type DMS ferromagnetism is governed by the kinematic mechanism involving the 
kinetic energy gain of heavy hole carriers caused by their hybridization with 3d electrons of Mn^^ 
impurities. Using the molecular field approximation the Curie temperatures Tc are calculated 
for several systems as functions of the impurity and hole concentrations. Comparison with the 
available experimental data shows a good agreement. 
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I. INTRODUCTION 



Dilute magnetic semiconductors (DMS) are semiconductors, in which transition or rare- 
earth metal atoms randomly replace a fraction of atoms in one sublattice. Transition metal 
(TM) impurities, due to their high abundance and diffusivity, can easily enter the host 
semiconductor. In the group IV elemental semiconductors, they occupy mainly interstitial 
positions, whereas in III-V semiconductors the 3d TM impurities usually substitute group 
III atoms. The metastability of the zinc blende phase of the DMS (III,Mn)V compounds 
and low solubility of manganese in these materials were the major obstacles for a synthesis 
of dilute magnetic semiconductors.^ However, the idea to combine in DMS the charge 
degrees of freedom of hole or electron carriers with the spin degrees of freedom of magnetic 
impurities became reality after a new doping technique based on non-equilibrium low 
temperature molecular beam epitaxy (LT-MBE)^'^ was developed. The discovery^ of 
ferromagnetic ordering in manganese doped InAs with the Tq = 7.5K Curie temperature 
fuelled DMS studies that resulted in fabrication of (Ga,Mn)As compounds with Tc = 
llOK^ik or even Tq = 140K.^ More recentlyi^Jl Curie temperatures exceeding 150K were 
reported in (Ga,Mn)As. An above room temperature ferromagnetism was announcediSiii 
also in GaN and p-tjpe GaP doped with Mn. Advanced III-V growth technique such as 
metal organic vapor phase epitaxy (MOVPE), or metal organic chemical vapor deposition 
(MOCVD) together with LT-MBE can produce good quality DMS with various element 
combinations. 

A unique coexistence of high-temperature magnetism and semiconductor properties 
opens new venues of fundamental studies and applications of DMS's. The DMS's may 
be quite promising in spintronics, information technology,^^*^^ quantum computing,— in 
manipulation of magnetism by an electric field^^ and/or by illumination,— in devices 
governed by the giant magnetoresistance effect,-^ where the interlayer tunneling resistance 
is changed under the action of a magnetic field. Among the variety of applications one 
cannot exclude other functionalities such as nanostructure aspect of DMS:— the DMS 
based quantum wells are expected to have unusual ferromagnetic properties with Tc (as 
shown in field effect transistorsi^) and high coercitivitiesi^ tunable electrically, rather than 
magnetically.-^ 

Although ferromagnetic order in (Ga,Mn)As has been observed for the first time as 
early as^l 1996 with Tq = 60K, its microscopic origin from the theoretical point of view 
still remains not well understood. Extrapolating from (Ga,Mn)As, FM ordering with high 
Tc was predicted for p-type (Ga,Mn)N2^ and found in Ref. 0. However, a direct ex- 
trapolation of the trends known for the GaAs-based material characterized by a relatively 
narrow forbidden energy gap to the wide gap GaN DMS is not well founded. The differ- 
ences in formation of the ferromagnetic order arise from the differences in the structure 
of chemical bonds between the Mn impurity and valence electrons in various III-V host 
semiconductors. The purpose of the present paper is to pinpoint these differences and to 
construct a general microscopic theory of the double exchange (superexchange) in dilute 
magnetic semiconductors. 

There are three approaches, which one may follow to study magnetic states of DMS's. 
According to the first approach, one chooses an effective spin Hamiltonian and calculates 
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the corresponding exchange couphng constant in the course of this derivation. The major 
part of the available descriptions of the FM order in such DMS's as (Ga,Mn)As, (Ga,Mn)P 
and (Ga,Mn)N, are based on semi-phenomenological models, which postulate the existence 
of local magnetic moments on the Mn sites, as well as an indirect exchange between these 
moments and the electrons in the valence band of the host semiconductor . i7i23^24^25^26 gQj^^g 
of these theories emphasize the role of shallow acceptor levels existing in the first two of 
these systems^! or of resonance levels.— 

The second approach is based on an extended cluster method, where a cubic supercell 
with a magnetic Mn ion in its center is used to calculate the density of spin polarized states 



of a "homogeneously doped" crystal (see e.g., Refs. l29ll3Clf ). The resulting picture provides 



information about positions and occupation of Mn-related majority and minority "bands" . 
Sometimes an effective s — d exchange Hamiltonian is introduced phenomenologically in 
this picture, and the effective exch ang e constant is chosen to fit the numerical data (see 
Ref. 2^ and the first paper in Ref. Isoi i. 

According to the third approach, realistic exchange parameters should be derived 
within the framework of a microscopic theory based on the knowledge of the electronic 
structure of III-V semiconductors doped with TM impurities. An exhaustive microscopic 
theory for isolated TM impurities in semiconductors was constructed more than two 
decades ago starting with the papers of Refs. 31 3^ (see monograph Ref. IH for a de- 
tailed description and Ref. |3J where a powerful numeric approach was developed). This 
theory contains all necessary ingredients for an accurate derivation of the indirect mag- 
netic exchange between interacting impurities. Our first attempt^ to apply this theory 
to (Ga,Mn)As resulted in a quantitative description of the hole concentration dependent 
Curie temperature Tc in p-tjpe DMS's within the framework of a very simple model 
Hamiltonian, in which the effective kinematic exchange results from hopping between the 
occupied d-states of the half-filled li-shell of the Mn impurity via empty heavy hole (hh) 
states near the top of the valence band. 

Here we present a general theory of the double exchange interaction in Mn-doped III- 
V semiconductors, which covers both p- and n-type materials. The theory establishes 
differences between the exchange mechanisms in these two cases and demonstrates that 
the conventional Zener double exchange mechanisn>2& is realized only in n-type DMS 
(presumably, in (Ga,Mn)N). 

It is well known that the TM ions may introduce impurity levels in the forbidden energy 
gap. Among them the manganese impurity in GaAs and some other semiconductors 
exhibit peculiar features. According to its position in the series of the TM elements, a 
neutral impurity Mn'^+ in the Ga-site is expected to have the 3d'^ configuration. However, 
the Mn ion retains its fifth electron in its 3d shell, because of an exceptional stability of 
the high-spin half-filled 3d^ state resulting from the strong intra-atomic Hund interaction. 
The manganese impurities not only introduce magnetic moments in the ^type (III,Mn)V 
compounds but also create potentials attracting holes, i.e. act as acceptors. Therefore, 
the neutral substitution impurity state is A'^{3d^p), where p stands for the loosely bound 
hole, and the manganese provides both holes and magnetic moments to the host, A is 
the spectroscopic notation indicating the irreducible re pre sentation describing this state 
(more about these notations can be found, e.g. in Ref. l33l ) 



3 



Existence of the complex Mn d^p in lightly doped bulk GaAs is detected by electron 
spin resonance measurement s2i and by infrared spectroscopy.™ These studies discovered 
an acceptor level inside the energy gap around 110 meV above the top of the valence 
band. This value is substantially higher than 30 meV resulting from conventional effective 
mass theory,— which already implies that treatment of the manganese impurities in III- 
V compounds requires a more refined approach. On the other hand, the presence of 
negatively ionized complexes A~{?)d^) is detected in (Ga,Mn)As epilayers grown by low 
temperature molecular beam epitaxy,^^ which indicates that both mobile and localized 
holes exist in ferromagnetic (Ga,Mn)As. The coexistence of weakly and strongly localized 
states is also in accordance with ac conductivity experiments.-^ 

Apparently, a similar impurity electronic structure is realized in (Ga,Mn)P. But in 
contrast, Mn in the wide gap GaN semiconductor releases one of its ci-electrons to the 
valence band, which is typical of all other TM impurities, and remains in the Mn3+(d'^) 
state, unless the sample is deliberately ra-doped. The latter state is neutral, so there is no 
need in binding a hole to maintain the neutrality of the ground state. This leaves no room 
for a direct extrapolation from (Ga,Mn)As to (Ga,Mn)N, and the kinematic inter-impurity 
exchange mechanisms in these two systems should be considered separately. 

In this paper we will develop a general theory of ferromagnetic exchange in DMS's. It 
turns out that the Zener mechanisn>2& in its classical form is realized in n-type (Ga,Mn)N, 
whereas the p-type systems (Ga,Mn)As, (Ga,Mn)P and, possibly, p-(Ga,Mn)N are exam- 
ples of another type of indirect (kinematic) exchange interaction, which was not observed 
in the family of Mn oxides described by the Zener theory of ferromagnetic insulators. The 
paper is organized as follows. Section II presents a detailed description of the microscopic 
model including the calculation of the exchange interaction between TM impurities. Since 
this calculation appears to be rather involved, its mathematical details are given in the 
Appendix. In Section III we apply the results obtained in Section II to ]o-type semiconduc- 
tors, find the kinematic exchange interaction and using the molecular field approximation 
calculate the Curie temperature as a function of the Mn content and hole concentration 
for j>-type (Ga,Mn)As and (Ga,Mn)P. The results are compared with available experimen- 
tal data. The n-type (Ga,Mn)N DMS's, which possesses specific features differing it from 
other DMS's, is discussed in Section IV. Finally, the scope of the theory and comparison 
with other approaches are summarized in the concluding Section V. 



II. INDIRECT EXCHANGE INTERACTION BETWEEN TWO MAGNETIC 
IMPURITIES 

According to the general principles of the theory of isolated TM in semiconductors ,??'^^ 
the electronic spectrum of these impurities in III-V semiconductors is predetermined by 
the structure of chemical bonds (hybridization) between the 3(i-orbitals of impurity ions 
and the p-orbitals of the valence band electrons, whereas their magnetic state is gov- 
erned by the Coulomb and exchange interactions within the Sd-atomic shell (intra-atomic 
electron-electron correlations), modified by hybridization with the electron states of the 
crystalline environment. 
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The localized TM impurity (i-states possess t2 and e symmetry of the crystal point 
symmetry group as revealed in the theoretical studies of electronic properties of isolated 



TM impurities in semiconductors carried out in the 70's and 80's (see, e.g. Refs. |41 
4^ 43..,44.45i ) . The e states are practically nonbonding ci-states of the impurity 3d shell, 



whereas the ^2 states form bonding and antibonding configurations with the p-states of the 
valence band. The two latter types of states are classified as the crystal field resonances 
(CFR), in which the d-component dominates, and the dangling bond hybrids (DBH) with 
predominantly p-type contribution of the valence band states. Crucially important is 
that the absolute positions of CFR levels relatively weakly depend on the band structure 
of the host semiconductor material. The CFR levels are pinned mainly to the vacuum 
energy, and this pinning is modulated by the counterbalanced interactions with valence 
and conduction band states.— The DBH states are more intimately connected with the 
peak in the density of states of the heavy hole band, which dominates the pd-hybridization. 
The variety of electronic and magnetic properties of Mn and other TM ions in different 
III-V host semiconductors stems from these universal chemical trends. 

Due to a large mismatch between positions of the valence bands in GaAs, GaP on 
the one side and in GaN on the other side, the Mn{d^ / d'^) CFR level falls deep in the 
valence band in GaAs, GaP and appears within the forbidden energy gap of GaN. The 
CFR nature of Mn [d^ /d'^) t2-level near the middle of the gap for (Ga,Mn)N is confirmed 
experimentally^ and follows also from numerical calculations.— This t2-level is empty 
in the neutral state yhx^^{d'^) of the substitution impurity. It becomes magnetically 
active only in n doped materials, e.g. in (Ga,Mn)N, where part of the Mn impurities 
capture donor electrons and transform into charged ions Mn^+(c?^). A similar behavior 
is characteristic of all other light TM ions in all III-V compounds including GaAs and 
GaP.^2^ 

Since the position of the valence band in GaAs and GaP is substantially higher in the 
absolute energy scale (with the vacuum energy as the reference point) than those in GaN, 
the CFR t2 level {d^/d^) in the (Ga,Mn)As and (Ga,Mn)P DMS's falls below the center 
of gravity of the valence band. Therefore, this CFR level is always occupied, Mn ions 
retain their fifth electron in the Sd-shell, and the neutral substitution impurity state is 
(3c?V), where p stands for a loosely bound hole occupying the relatively shallow acceptor 
level in the forbidden energy gap of GaAs and GaP semiconductors. 

Both of the above mentioned situations may be described within the resonance scat- 
tering model of the TM impurity states in semiconductors based on the Anderson single 
impurity Hamiltonian, proposed originally for TM impurities in metals,— and later mod- 
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ified for semiconductors in Refs. 31 3 An extension of the Anderson model to 



two TM impurities in metals was proposed in Ref. |42l resulting in a ferromagnetic cou- 
pling in TM doped metals. Using this model Ref derived an interaction similar to the 
RKKY interaction. Here we discuss a model describing interaction between Mn (and other 
TM) impurities in a semiconductor host. In accordance with Hund rule for the electron 
occupation in Mn ions in a tetrahedral environment,— >^^^ two competing states Mn(d^) 
and Mn (d^), involved in the double exchange in DMS, have the d'^{e^t^) and d^{e^t^) 
configurations, respectively. The next charged state d^{e^t^) has a much higher energy 
due to the strong intra-ionic Coulomb interaction. Therefore, the 'passive' nonbonding 
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e-electrons may be excluded from our consideration. The indirect exchange between mag- 
netic moments is mediated by virtual transitions of ^2 electrons into unoccupied valence 
band states. 

The Hamiltonian for two TM impurities in a III-V semiconductor has the form: 

H = Hf, + H, + Hm, (1) 

where the band Hamiltonian 

HH = J24ik.c^^'^ (2) 

p.cr 

includes only the heavy hole band. Here c^phai^pha) is the creation (annihilation) operator 
of a hole with momentum p and spin orientation a in the hh band of the semiconductor 
with the energy dispersion e^. The heavy hole band gives the dominant contribution to 
the formation of the impurity states , ^^i^'^ and governs the onset of ferromagnetic order. 
The second term in the Hamiltonian describes the electrons within the Mn atoms 
with a possible account of the crystal field of the surrounding atoms. In principle, we 
have to write down a multielectron Hamiltonian describing the degenerate states of the 
d shell, and including Coulomb and exchange interactions. However, as we show below it 
is usually sufficient to consider the non-degenerate version of this Hamiltonian 

Hd = ^ (Edhia- + ^hiafii.oA , (3) 

which simplifies the calculation considerably. Here is the atomic energy level of the lo- 
calized Mn t2-electrons and \J is the Anderson-Hubbard repulsive parameter; rijo- = dl^dicj 
is the occupation operator for the manganese t2 - electrons on the impurities, labelled 
i = 1, 2 in (III,Mn) V DMS. The situations when the degeneracy of the ^2 states is impor- 
tant and the exchange interaction (Hund rule) in the multielectron atom starts playing an 
essential part will be outlined and the appropriate corrections will be introduced, when 
necessary. 

The last term in Eq. describes the scattering of heavy holes by the impurities 

Hhd = Hj[^ + HjfJ , (4) 



p. o-j' 

^hd — 2-^ ^^PP'^pha^P'ho-fi 



Here Vpd is the p-d hybridization matrix element, and Wppi is the scattering matrix 
element due to the difference of the pseudopotentials of the host and the substituted 
atoms. The direct overlaps between the d-electron wave functions of the neighboring 
impurities in the DMS and the Coulomb interaction between them are neglected. 

As mentioned above, the orbital degeneracy of the t2 states is not taken into account 
in this simplified model. In reality, the orbital quantum numbers are important at least 
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in three respects. First, the half-filled subshell of Mn is occupied by e and t2 electrons 
with parallel spins in accordance with the Hund rule, so the sixth electron can be captured 
in the configuration only with the opposite spin. The energy cost of this capture is 

U + J, where J <C f/ is the exchange energy. This feature of the impurity is preserved in 
the above simplified Hamiltonian: Q the reaction d^ + e —>■ d^ is changed for d^ + e —>■ d"^, 
with the spin of the second electron opposite to the first one. The energy cost of this 
reaction is U, and the principal features of the ion, important for the formation of the 
localized magnetic moment, are practically the same as in the original atom. 

Second, it is the Hund rule that requires that the total angular moments of the Mn 
atoms be parallel to allow the indirect inter-impurity interaction between the high-spin 
states via the hh valence band states. We will take the Hund rule explicitly into account, 
when calculating Tq, and it will be shown below, that the energy gain due to this indirect 
interaction, monitored by the Hund rule, leads eventually to ferromagnetic ordering in 
DMS's. 

Third, the three-fold degeneracy of the ^2 electrons is also manifested in the statistics 
of the localized states, and therefore it infiuences the position of the Fermi energy in the 
energy gap. These statistics are described in the Appendix. The degeneracy introduces 
also numerical factors in the effective exchange constants. These factors are also calculated 
in the Appendix. 

Having all this in mind, we proceed with the derivation of the indirect inter-impurity ex- 
change within the framework of the non-degenerate Alexander-Anderson model, whereas 
the corresponding corrections for the real orbital structure of the Mn ion will be made 
when necessary. To calculate the energy of the exchange interaction between two magnetic 
impurities, one should find the electronic spectrum of the semiconductor in the presence 
of two impurities. Each impurity perturbs the host electron spectrum within a radius 
Th. Inter-impurity interaction arises, provided the distance between the impurities Rij is 
comparable with 2rfe. General equations for the two-impurity states and the correspond- 
ing contribution to the exchange energy are derived in Appendix. Here we present only 
the final equations, which will be used in the derivation of the effective exchange coupling 
and the Curie temperatures for the specific DMS's. 

The quantity, which we actually need here is the impurity related correction to the 
energy of the system. It is given by the standard formula (see, e.g., Ref. 00 

1 r°° 

^magn ^ / ^ . TrAG[£ - iSsign{e - fx)]de - ^ Undi^Udn. (5) 

AG = G(£:) — g°(e) is the two impurity correction to the total one-particle Green function 
G(e) = (e-H)-^ 

Here G is defined as a matrix in {h, d) space, and g° is the same matrix in the absence of 
impurity scattering, see Eq. ()A8|) . The Green function is calculated in the self-consistent 
Hartree-Fock approximation for the inter-impurity electron-electron repulsion U, which is 
sufficient for a description of magnetic correlations in TM impurities in semiconductors. 
This solution is described in the Appendix. Then the part of the two-impurity energy ^ 
responsible for the inter-i mpurity magnetic interaction may be found (see Eqs. ()A17|) and 
()A18|) ). As shown in Refs. I3lll45l the resonance scattering alone may result in creation of 
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CFR and DBH levels split off from the valence band, provided the scattering potential is 
comparable with the bandwidth. However, in the case of shallow DBH states the potential 
scattering may dominate their formation. Like in the well-known Koster-Slater impurity 
model, the potential scattering in our model is described by a short-range momentum 
independent coupling constant, Wpp' ~ W, and the same approximation is adopted for 
the resonance scattering parameter, Vpd ~ V. 
Then (0) acquires the compact form 



Here the first term in the r.h.s. describes the band contribution to the exchange energy, 
and the integration is carried out from the bottom of the hh band {£hb) to the Fermi level 
Ep- In this integral the kinetic energy gain of the band electrons due to scattering by 
ferromagnetically aligned impurities is incorporated. The contributions of resonance and 
potential scatterings are given by the first and second term in this integral, respectively. 
The last term AEioc is the contribution of the occupied localized CFR and/or DBH states. 
These states are described by zeros of the functions R'^(e) and Q'^(e) in the discrete part 
of the energy spectrum (see Eqs. ()A11|) and ()A9|) . respectively). 

In all the cases the energy gain results from the indirect spin-dependent inter-impurity 
overlap, and the mechanism of the effective exchange interaction may be qualified as kine- 
matic exchange. In general, double exchange favors a FM order, because the splitting of 
the energy levels belonging to two adjacent impurities occurs due to electron hopping 
via unoccupied band or impurity-related states without spin reversal. Were the impurity 
angular moments non-parallel, the Hund rule would suppress the probability of the elec- 
tron, with spin parallel to the first impurity angular momentum, to hop onto the second 
impurity with a "wrong" direction of the angular momentum. 

The level splitting, when not suppressed by the Hund rule, results in an energy decrease 
provided that not all of the available band and impurity levels are occupied. One should 
note that this kinematic exchange cannot be reduced to any conventional double exchange 
mechanism because of an actual interplay between three contributions to the magnetic 
energy, namely, the scattered valence band electrons, the CFR states and the DBH states. 
We start with the analysis of the band contribution AEb^ex- Since in all the cases we 
deal with nearly filled hh bands, it is convenient to calculate the energy in the hole 
representation. 

By means of a simple trick (as described in the Appendix) the contribution of the hh 
band in the basic equation (jHl) may be transformed into 



where Eht stands for the top of the valence band. This equation is obtained from the more 
general equation ()A21|) by neglecting the potential scattering in the band continuum. 
Then, inserting Eq. ()A11|) for the matrix R'^ie) we obtain 




(6) 




(7) 



AE, 



'b,ex 




L 



de In 



gf{e-i5)-V^L^i{e-i5) 
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Here 

Li.{e-t6) = Pi,j{e) + '-rij{e) (9) 

are the standard lattice Green functions ()A10|) of the continuous argument e. Retaining 
only the leading (quadratic) terms in Lu in Eq. (jHl), the exchange energy due to the hh 
band becomes 

A^,.. ^ r ,e ^i(£)lH(£) . (10) 

' \e-E,-V^P,,(e)Y + —T\,(e) 

Eq. (jH)) with Li2{e) 7^ holds only if the impurity angular moments are parallel to 
each other (ferromagnetic case), since the electron hopping between the impurities via 
the hh band takes place without spin-flips. If the angular moments are non-parallel, 
the Hund energy does not allow the electron to hop to the neighboring impurity with 
a "wrong" angular momentum direction. So the impurities with non-parallel angular 
moments cannot exchange their electrons, which effectively means that for them L12 = 
L21 = 0. Thus, Eq. (fTn|l is the energy decrease due to the kinematic indirect exchange 
between a pair of magnetic impurities with parallel angular moments. If the moments are 
not parallel the corresponding energy decrease is zero. 

Another contribution to the magnetic energy AE (Eq. (jHl)) originates from the inter- 
impurity hopping via empty localized states, if available. As was mentioned above, there 
are two types of such states, namely CFR and DBH-type levels. For the CFR states, one 
may, to a good approximation, neglect the potential impurity scattering. To take properly 
into account the Coulomb blockade effect on the impurity site, one has to calculate the 
Green function in a self-consistent way, known as "Hubbard I" approximation (see Eqs. 
flMIll . (|XT6|l in the Appendix). 

As was pointed out above, the origin of the energy gain is the inter-impurity level 
splitting. The latter is derived in the Appendix (Eqs. ()A22|1 and ()A23|l ). The result is 

,E.,,,^±A,'JL_I_:^, (11) 
where the splitting of the levels of the isolated impurities is 

The second term in Eq. results from the repulsion of the two-impurity levels from the 
band continuum. All the functions Pij and their derivatives P[^ are taken at e = E^p^, 
i.e. at their positions for the isolated impurities ()A22|1 . 

Considering the contribution of the DBH states we may apply the procedure outlined 
in the Appendix and to find the kinematic exchange due to the DBH states lying above 
the top of the valence band 

KV'^P'^ 
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Now all the functions in Eq. ()13p should be calculated at the DBH level {e = E^^^). This 
part of the kinematic exchange favors the ferromagnetic ordering, since A'i(-^d_b//) < 0- 
If the DBH levels lie not too far from the top of the valence band their contribution may 
be comparable with that of the hole pockets and should be properly taken into account 
when calculating the Curie temperature. At a sufficiently high Mn content the splitting 
of the DBH level may result in the formation of an impurity band and its merging with 
the hh band. This case is discussed in the next section. 

Thus one sees that the contribution, AEioc, of the localized states to the kinematic 
exchange is not universal. It depends on the type of conductivity and should be considered 
separately for p- and n-type materials. In the p-type samples one should take into account 
the empty shallow Mn-related levels, which are present in p-(Ga,Mn)As and p-(Ga,Mn)P 
and apparently absent in p-(Ga,Mn)N. In the n-type (Ga,Mn)N, the localized states are 
due to the deep Mn CFR levels, and hopping over these levels alone determines the 
exchange mechanism. The double exchange mechanisms in these two cases are obviously 
different. 

To calculate the Curie temperature Tc for the p- and n-type materials one should 
extend the two-impurity calculations for dilute materials with small Mn concentration x 
and take into account the corresponding transformations in the magnetic energy. This 
program is realized in the two next sections. 

III. MAGNETIC ORDER IN P-TYPE DMS 

Here we discuss the formation of ferromagnetic order in Mn-doped III-V p-type semi- 
conductors. (Ga,Mn)As together with (Ga,Mn)P are the most celebrated among them. 
As mentioned above, the CFR level of the fifth d electron lies in these systems below 
the heavy hole band. A DBH level is also formed above the top of the valence band. 
At sufficiently high concentrations these DBH levels start broadening into an impurity 
band, and may merge with the valence band. The band structure of (Ga,Mn)As is shown 
schematically in Fig. ^ 

The empty localized levels lying above the Fermi energy. (While making fitting of 
the theory to experimental data we use the Fermi energy ep instead of the chemical 
potential /i, since they nearly coincide under the condition of the experiment.) appear 
due to the combined action of both the potential (W) and resonance (V) scattering 
mechanisms. Formally they can be found as zeros of the determinant R(£:) (see Eq. (jAlljl ). 
Neglecting small corrections due to resonance scattering, the energy of an isolated DBH 
level, E^j^^fj = Ei, corresponds to a zero of the function 



like in the Koster-Slater one-impurity problem (see the graphic solution on the left panel 
of Fig. P). The calculations below will be carried out assuming the model semi-elliptic 
density of states for the hh band 



q{ei) = 1 - WPii{ei) = 0, 



(14) 
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FIG. 1: Left panel: graphical solution of Eq. (jlTj) for the bonding CFR and Eq. ()14() for 
the antibonding DBH levels. Right panel: energy levels in (Ga,Mn)As. The CFR d-levels 
d^/d^ (denoted by Ri(2)) of each impurity, lie below the hh band. The DBH levels (energies eio-, 
£20-) are split from the hh band and form localized (acceptor) levels in the energy gap. The CFR 
levels d^/d^ R_i(_2) lie high in the conduction band. The impurity band is shown by the shaded 
area together with the position of the Fermi level as a function of the width of the impurity 
band (horizontal axis) 

where w is the hh bandwidth and the energy e is counted from the top of the hh band. 
Then the position of the acceptor DBH level can be found explicitly 

^. = Kb„ = W (1 - jj^) . 

Using a realistic value w = 2.9eV— for the hh bandwidth and taking the matrix element of 
the potential scattering W = 1.02 eV, one obtains the acceptor level position Si = 85meV 
for (Ga,Mn)As. The influence of resonance scattering on the positions of the shallow DBH 
levels of the two interacting impurities may be usually neglected (see Ref . for a detailed 
description of the interplay between resonance and potential impurity scattering). 

At a sufficiently high impurity content, x, an impurity band is formed in (Ga,Mn)As 
since the DBH levels are split due to the indirect interaction within the impurity pairs. 
We neglect the pd-hybridization V, keep only the contribution of the potential scattering 
(W) and obtain 

Se^ = -^q±, (15) 
where q_|_ are the two solutions of Eq. ()A25|) . 



11 



Then Eqs. fll5|) and ()A25|) result in 



2' ^ Pn' 
The functions P and P' are calculated at e = ^j. The half- width of the impurity band 
can be roughly estimated as zAdbh{R) where R is the typical distance between the 
impurities, which depends on the impurity content x. The z value characterizes the 
number of neighboring impurities participating in the interaction with any given impurity. 
If the half- width zA£)bh{x) of the impurity band exceeds the energy Si (counted from 
the top of the valence band), the impurity band merges with the valence band and they 
both form a unified continuum of states. For the above mentioned values of the model 
parameters this merging occurs even at x < 0.01. 

Eventually, the magnetic order is due to the exchange interaction between the occupied 
CFR levels. These levels correspond to the states d^/d^ of the Mn ions below the hh 
band, whereas the empty d^ /d^ CFR levels are pushed up to the conduction band by the 
Anderson-Hubbard repulsion U (Fig. H]). 

The left panel of this figure depicts a graphical solution of the equation 

R{e) = 0, (16) 

with R(e) defined in Eq. ()A11|) . Neglecting the interaction between the impurities, L12 = 
(and, hence Mi2{R) = 0), Eq. (fTH|) takes the form 

EcFR = Eia- = Ed + V'^Lii{Eio-), (17) 

and describes the formation of the deep impurity states of CFR type out of the atomic 
d-levels below the hh valence band. The values of the parameters used in the graphical 
solution presented in Fig. [T]will be discussed below. 

For a finite distance between the impurities, the overlap of impurity wave functions 
due to the nonzero Li2{R) leads to a level splitting as shown in the right panel of Fig. 
Each of the two impurity levels is singly occupied due to the Anderson-Hubbard on-site 
repulsion. The inter-impurity overlap arises due to the pd-hybridization, which does not 
involve spin-flips, hence it occurs only for the parallel impurity angular moments (Hund 
rule). There is no level splitting for the non-parallel impurity angular moments, unless 
one takes into account indirect interaction via empty d^ states. 

FM ordering arises, provided the state with the parallel impurity angular moments is 
energetically preferable in comparison to that with the non-parallel ones. It is obvious that 
the splitting per se cannot give an energy gain, when both states are occupied. One should 
take into account all the changes in the energy spectrum, namely, the reconstruction of 
the partially filled merged impurity and hh band. The indirect interaction involving empty 
states near the top of the hh band is in fact a novel type of the double FM exchange, which 
resembles the well known Zener exchange^ but differs from it in many important aspects 
(see below). 

The contribution favoring the FM order can be obtained from Eq. ()10|) with the 
addition of the part due to the impurity band merged with the valence band, 

^E^M = -- de ^M^)PMe) 

^ [e-E,- V^Puie)? + ^Ff, (5) ' 
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(the reference energy in this equation is taken at the top of the hh band). In the FM- 
ahgned spin configuration only the majority spin subband contributes to AEpM, therefore 
the spin index is omitted for the sake of brevity. 

We use in our numerical estimates the following approximate relations, 



e — u J k (iu) R e — u 

Ti2{e) ~ 2Txp{e) sm[k {e) R]/k {e) R, Tnis) ^ np{e) = T22ie). (19) 

(here R is the inter-impurity distance). The value of the wave- vector k is found from the 
equation e = ehh{k), where Shhik) is the hh energy dispersion. In our model calculations 
the impurity band merges with the top of the valence band of (Ga,Mn)As at the concen- 
tration Xcrit = 0.0065 with the renormalized 3d-wave-function enveloping 5 to 6 unit cells. 
Earlier transport measurements for low Tc samples indicate merging even at x = 0.002. 
All this justifies the approximation adopted for the hole state in Eq. (fTHjl . This equation 
is our working formula, from which we obtain Tc- 

The mechanism outlined above competes with the antiferromagnetic (AFM) Anderson- 
type superexchange-ft involving the empty (f' states. The energy gain in the latter case is 
estimated as 

AEafm ^ -V^LlJU . (20) 
It should be compared with the above Zener-type mechanisms. We assume that 

Pij ~ w-\ Ti, ~ 4^'u;-3/^ ep-Ed- V^Pniep) = ^aV^w , 

with a < 1 (see left panel of Fig. Q). Then one finds from Eqs. (fTHj) and 
((20)) that \AEafm\ ~ V^a'^exp{-2KhRi2) /{Uw'^Rl^), with Kh = ^2m{eb^ - Eda)IK 
\AEfm\ ~ 2eF {sf /w)^^'^ / \{^Cif' + I {^w)\ and FM coupling is realized provided 
\AEpm\ > \AEafm\, which normally takes place since the parameter U is large. Nu- 
merical estimates are given below. 

Now the difference between the double exchange mechanism in DMS and the Zener 
double exchange in transition metal oxides becomes clear. The conventional Zener double 
exchange mechanism was proposed^^ for (La, A^"'")Mn03. In this case Mn ions are in 
different valence states (Mn^"*" and Mn^+). In our case it would have meant that one 
of the two levels Ecfr] were empty. It is seen from Fig. [T] that in (Ga,Mn)As both 
these states are occupied, and the Zener mechanism in its original form is not applicable. 
Actually, the FM order in the p-tjpe DMS arises due to the kinematic exchange between 
the two Mn(d^) ions via the empty states near the top of the valence band. 

In order to calculate the exchange energy Eq. ()18|) . the dependence of the Fermi level 
on the Mn content eF{x) should be known. It is found from the equation 

5w{x) 

Xs = 2 [ p{e)de, (21) 



Sw{x) = Ei + zAf)bh{x), and the hh density of states is approximated as 

p (e) = 8/ [vr {w + 6w{x)f] ^/[6w{x) -e]{e + w)e [5w{x) -£\9{£ + w). 
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{z is the coordination number in the impurity band). The per site hole concentration Xg is 
proportional to the hole density ph] Xg = a^ph/8 in the zinc-blende structure. Equations 
(|18p and (|2H) allow one to determine the pair exchange energy as a function of Mn con- 
centration X and connect it with the experimental data for the given hole concentrations 
Ph{x) which are taken from the measurements.— 

The Curie temperature Tq was calculated in the molecular field approximation. Ac- 
cording to this approach the spin Hamiltonian reads 

Hmf = 2 ^EFM(Rij)Ji ■ (J) , (22) 

i 

where summation runs over all positions of the Mn impurities with an angular moment 
operator Jj. Here the factor 1/2 accounts for the fact that the energy gain of the FM vs 
AFM orientation of two coupled spins is AEpMCR-ij), contrary to 2JHeis in the Heisenberg 
model. Then the Curie temperature can be readily found, 

AEpMjx) zJ{J+l) 
ks 6 

where 1 is, similarly to z, a measure of the neighboring atoms participating in the exchange 
interaction. It is certainly close to z, although it does not necessarily coincide with it. It 
is worth mentioning here that unlike the case of magnetically doped metals^ we cannot 
represent the effective inter-impurity coupling in the asymptotic form 

AEfm{R^j) ~ cos{2kFRij + (t))/{kFRijf 

because typically /c^^ ~ FRij in doped semiconductors. 

In (Ga,Mn)As the total angular moment of a complex Mn(3(i^p) is unity^S: J = 1, 
since it is formed by the moment j = 3/2 of the loosely bound hole AFM-coupled to the 
Mn center with the spin 5* = 5/2. Then the numerical factor in Eq. is close to unity. 
The results of our calculations are presented at Figs. |21 and El 

In order to find the dependence Tc{x), the calculated value of the Fermi level via 
Eq. (jSH), based on the experimental data of Ref. EE is used in Eqs. p3|l and (|T8|l . A 
polynomial fit is used for the hole density ph{x) in Eq. (|TH|l . The values of the model 
parameters characterizing the impurity d-state are U ~ 4.5 eV, V = 1.22eV, while the 
hh mass m = 0.51 • rriQ and hh bandwidth w = 2.9eV were taken from Ref. IH^. The 
hybridization strength V was obtained from Eq. (fT7|l with Ecfr = — 3.47eV {E^pj^ = 
— 3.4eV according to Ref. IsgI ). The value of Ei = 85 meV was chosen for the energy of 
the shallow acceptor level (not far from the experimental value^"^ ef^^ = 110 meV.) The 
value of z = 2.5 was taken for the coordination number in as-grown samples (Fig. 2, lower 
curve). At these values for the model parameters the ratio \ AEfm\/\A.Eafm\ ~ 2 justifies 
the dominance of the FM coupling in (Ga,Mn)As. 

The non-monotonous dependence Tc{x) is due to a non-equilibrium character of the 
sample preparation. Apparently, the ratio between the concentration of Mn impurities and 
the actual hole concentration depends on the doping method and the thermal treatment. 
In particular, the annealing of the sample results in a reduction of the donor-like Mn- 
related interstitial defects in favor of acceptor-like substitution impurities . To describe 
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FIG. 2: Calculated dependence of Tc on the manganese concentration x in (Ga,Mn)As based on 
the experimental data for the hole density^^ ph (x). Solid squares (circles) stand for experimental 
Tc {x) of annealed (as-grown) samples. Broken lines take into account the error bars of the hole 
density. 

the annealing effect we changed the value of z from 2.5 to 4 in Eqs. p8|) and ()23p . The 
results of numerical fitting are shown in the upper curve of Fig. |21 

Recent detailed measurements of the hole concentrations in a series of both as-grown 
and annealed Gai_a;Mna;As samples,— allowed us to compare the theoretical plot Tc{p) 
with the experimental data. These results are presented in Fig. El Our fitting procedure 
uses the same equations and (fTS|). and the same values for the model parameters 
EcFRi U, Si as in the above estimate. The Mn concentration point x = 0.067 is used as 
a reference point, and the coordination number ^ = 4 is chosen. Two theoretical curves 
correspond to two values of the hybridization parameter V = 1.22 and 1.24 eV (solid 
and dash-dotted curves, respectively). One may conclude from these two fittings that 
the theory is not very sensitive to the choice of the model parameters. To check this 
statement, we made one more fitting of the experimental data obtained only for annealed 
samples, (see Fig. These data are taken from Ref. IH^ We see that the calculations 
with the same set of model parameters give satisfactory quantitative agreement with the 
experiment for these samples as well. 

Thus, our theory provides a satisfactory quantitative description of the behavior of the 
Curie temperature Tq as a function of the Mn content x and the hole concentration 
in p-(Ga,Mn)As. This description is applicable also to p-(Ga,Mn)P. The band structure 
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FIG. 3: The dependence of the Curie temperature on the hole density p in (Ga,Mn)As. Sohd 
and dash-dotted Unes are the theoretical curves for the two values of the hybridization parameter 
(V = 1.22 and 1.24 eV, respectively). Experimental points (filled squares) are taken from Ref. 
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of (Ga,Mn)P can be also schematically represented by Fig.l, however, the DBH level lies 
deeper in the energy gap than in p-(Ga,Mn) As: its position is estimated as ~ 0.4 eV above 
the top of the valence band (see, e.g. Refs. .33m34m47h52ii) . The impurity band appears to 
be narrower than in (Ga,Mn)As. It does not merge with the hh hand at any reasonable 
Mn concentration and its contribution to the indirect exchange can be neglected. We 
tested our theory by fitting the experimental valued of Tq = 270 K reported for a p-tjpe 
(Ga,Mn)P at x = 0.06. The impurity level Si = 400 meV arises at W = 1.4 eV. We 
have found that the value of Tc = 270K is reproduced by means of Eq. ()23|) with the 
following fitting parameters: w = 2.6 eV, V = 1.31 eV, x = 0.05, z = 4, E^i = —1.3 eV 
(then EcFR = —3.28 eV). The chosen value of ph = 1 ■ 10'^^ cm~'^ is the hole density in a 
C-doped GaP before Mn implantation.^^ So the values of our model parameters are close 



to those for (Ga,Mn)P in the experiments of Ref. |5J. Besides, these parameters are in 
a good correlation with the chemical trends in the systems (Ga,Mn)As, (Ga,Mn)P (the 
latter has a smaller lattice spacing and larger effective mass m^^^ than the former). 

As for |?-(Ga,Mn)N, the theory should be somewhat modified in order to describe 
this material. In this case Mn remains a neutral substitution isoelectronic defect in the 
configuration Mn^+((i^), since the Mn^+/^+ transition energy falls deep into the wide 
energy gap in accordance with the experimental^ and theoretical^ data. If the hole 
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FIG. 4: Calculated Tc (triangles connected by a solid line) vs the annealing temperature Ta 
in (Ga,Mn)As. The data for the hole density and the experimental values of Tc(Ta) (closed 



squares) are taken from Ref. 
z = A.O,EcFR = -2.6eV. 



58l . The parameters used in the calculations are: V = 1.22eV, 



concentration, p^, exceeds the Mn content, x, then we return to the situation discussed 
in this section, but with the locahzed moment J = 2 characteristic of the Mn^+(e^t^) 
configuration. In the opposite case, the valence band is full, and the Fermi energy lies 
within the deep impurity band formed by the Mn^+Z^"*" levels. This situation is discussed 
in the next section. 

Unfortunately, there are not enough experimental data on Tc{x) available to make a 
detailed comparison with the theoretical predictions both in p-(Ga,Mn)P and j»-(Ga,Mn)N 
materials. Moreover, no p-type conductivity was observed in (Ga,Mn)N even in the cases 
when the pristine GaN crystals were of a p-type. Due to the lack of experimental informa- 
tion, we have confined ourselves with a quantitative description of the p -type (Ga,Mn)As 
and turn now to the case of n-(Ga,Mn)N, where the double exchange mechanism, as 
described above, should be revisited. 



IV. MAGNETIC ORDER IN iV-TYPE DMS 

In accordance with the general predictions of the theory , Mn is a deep acceptor in 
GaN, and each Mn impurity creates an empty t2 level close to the middle of the energy 
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gap of GaN. Therefore, the Fermi energy should be pinned to these levels unless other 
(shallow) acceptors create enough free holes near the top of the valence band. This 
statement is confirmed by extended cluster "quasi band" calculations.— If a Mn-doped 
sample contains also shallow donors in a noticeable concentration, then the deep Mn- 
related impurity band is partially filled by electrons, and one arrives at the problem of 
magnetic order in n-type (Ga,Mn)N. 

In this case the d^/d!^ CFR energy levels forming the impurity band are partially 
filled, so one encounters the mixed valence Mn^^/Mn^^ situation similarly to the original 
Zener model.— Hopping in the impurity band is possible due to an overlap of the tails of 
the impurity wave functions. These tails are formed by the superpositions of the Bloch 
electron wave functions from the hh band, so the latter in fact plays the role which is 
similar to that of the oxygen p-orbitals in La(Mn,Sr)03 considered by Zener. 

The specific feature of the n-type (Ga,Mn)N DMS's is that the d^/d'^ CFR impurity 
levels, EcFR-i lie deep within the forbidden energy gap. At a high enough Mn content they 
may form an impurity band, which still remains well separated from the hh band. This 
impurity band may be partially occupied by electrons and be responsible for the possible 
magnetization of the DMS. To describe formation of a deep impurity band within the 
framework of our model, we solve Eq. ()17|1 for a Mn-related CFR level and then substitute 
this solution into the secular equation (fT^ . These calculations are performed by using the 
same semi-elliptic density of valence states p{e) as in Section III. The bare level Ed in this 
case is above the top of the hh band, and the pd hybridization only slightly renormalizes 
its position in the energy gap. The graphical solution of Eq. (fTTjl is presented in the left 
panel of Fig. 

We assumed that the hybridization parameter is nearly the same as in (Ga,Mn)As, 
V = 1.2 eV. The right panel of Fig. El illustrates the formation of an impurity dimer 
via an overlap of the tails of the CFR wave functions in accordance with Eq. (fTT)|) . 
Then the electron hopping between the correlated states -Ri and R2 initiates formation of 
an impurity band and the Zener type indirect exchange interaction arises, provided the 
impurity band is partially occupied. In principle, the double exchange via empty states 
in conduction band should be taken into account. However, this contribution is small 
in zinc blende host crystals, because the hybridization Vds between the impurity (i-states 
and the s-states of the electrons near the bottom of conduction band is weak due to the 
symmetry selection rules.— The empty (d^/d^) CFR levels are pushed up high to the 
conduction band by the Anderson-Hubbard repulsion U. Like in the previous case these 
levels can be a source of AFM ordering in the impurity band, but we do not consider this 
relatively weak mechanism here. Another contribution to indirect exchange in the case of 
n-(Ga,Mn)N is the fourth-order in |Vrfpp|Vdsp Bloembergen-Rowland indirect exchange.-^ 
However, this term involves the hybridization with both valence and conduction electron 
states, so it is small due to the above symmetry ban. 

In order to calculate the double exchange interaction in such an impurity band, we 
consider first a pair of impurities at a distance R from each other. To be definite we 
assume that the chemical potential fi lies below the level Ecfr of an isolated impurity. 
According to Eq. ()12|) the interaction between the impurities results in splitting of the 
impurity levels with half-width A(i?), where we neglect its possible dependence on the 
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FIG. 5: Left panel: graphical solution of Eq. ((T7jl for the antibonding CFR level in (Ga,Mn)N. 
Right panel: energy levels in Ga(Mn)N. The CFR d-levels d^/(f (denoted by Ri(2)) of each 
impurity, lie within the energy gap. The CFR levels S /d^ R_i(_2) lie high in the conduction 
band. 

direction of the vector R. 

It is emphasized here again that the pd hybridization does not involve spin flips. It 
means that the electron spin does not change its projection in the course of the whole 
process leading to the indirect exchange, during which the electron is transferred from 
one impurity to the other and back via band states. This process is possible only if 
the spins S, of the two impurities are parallel. Otherwise the Hund rule will require an 
additional energy when the electron tries to hop onto the impurity atom with a "wrong" 
spin orientation. 

One may conclude that the indirect exchange interaction within a pair of impurities 
appears if their angular moments are parallel (FM ordering), and if the distance between 
the impurities is not too large, i.e. 

R < i?ex(/i) 

where Rex{lj) is the solution of the equation 

A(i?e.)=i?SF«-/^- (24) 
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Then one of the two impurity levels sinks below the chemical potential /i and is occupied 
by an electron resulting in the energy decrease A{R), whereas the second level remains 
empty and does not contribute to the energy balance. If the angular moments of the 
impurities are not parallel, Hund rule forbids the indirect exchange and there is no energy 
decrease. Therefore, we conclude that the exchange energy in such a pair equals just 
—A[R). It is zero if at least one of the above two conditions is not met. 
The pair exchange interaction can be estimated as 

A = / A{R)g{R)dR (25) 



^0 

where g{R) is the impurity pair distribution function. For example, we may assume that 
the random distribution of the impurities obeys the Poisson law 

g{R) = ^exp\-^A (26) 



3 ^ [ R\ 

where R is the average distance between the impurities. The latter is connected with the 
impurity content x by ^R = d^, where do is the minimal possible distance between two 
Mn ions. This assumption may not work too well for relatively high impurity content, since 
it overestimates the contribution of closely lying impurities. In fact the impurities cannot 
approach each other to distances smaller than the distance between two neighboring Ga 
atoms, a-\/2/2, a lattice constant. In such a case corrections should be introduced for 
these small distances. 

In order to complete our calculations, we need also an equation connecting the position 
of the chemical potential with the electron concentration n in the impurity band. This 
concentration coincides with the concentration of impurities in the d^ state, hence 

n= p{e)nds{e — fi)de (27) 

where p{e) is the density of states in the impurity band. The distribution n^s {e — /i) is 
calculated in the Appendix, Eq. ()A16|) . 

Since, in unfilled impurity band the chemical potential and the Fermi level ep practi- 
cally coincide, we use below the latter term as in the case of p-type samples instead of 
/i. At low temperatures nj_5{e — ep) is 1 for £ < ep and for £ > ep. This allows us to 
circumvent the calculation of the density of states p{e) and present the condition ^I7\\ in 
the alternative form 

n= g{R)dR, (28) 



^0 

which is more convenient for our calculations. If an experiment provide us with depen- 
dencies of the Curie temperature and the electron concentration on the Mn content, then 
Eqs. fl24|) and ()28p will allow us to connect the Fermi energy with the electron density in 
the impurity band and after that to use Eq. in order to calculate the average pair 
exchange interaction A(x) as a function of the impurity content x. 

We may now use the theory of dilute ferromagnetic alloys developed in Ref. [gO (see 
also the review Ref. |6l|) and arrive at the final result of this procedure, i.e. the Curie 
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temperature. Unlike the case of (Ga,Mn)As, the magnetic energy in impurity band is 
determined by the spin S of Mn 3c? shell, which equals 5/2 and 2 for and c?^ config- 
urations, respectively. In our crude approximation it is sufficient only to know that the 
proportionality coefficient between Tc and A is of order of unity: 

6kB 

At low electron concentrations, when Rexisp) > do-, only a minor part of the impurities 
is coupled by the indirect exchange interaction. Using the above averaging procedure 
for calculating Tc does not hold any more and a percolation type of approach should be 
applied. The resulting Curie temperature may become low and will decrease exponentially 
with the concentration.i^ 

The behavior of the Curie temperature for higher electron concentrations, when the 
Fermi energy lies above the energy -E^'fr' "^^^ be found in a similar fashion. We just have 
to switch to the hole representation and use the distribution function ^iridi^e — Sf) instead 
of niff>{e — Sp)- (Dealing with holes we have now to account for the three-fold degeneracy 
of ^2 states.) The function 3n^4 is 1 ii e > Sp and Q ii e < Sp- As a result, we will obtain 
a mirror symmetric dependence of Tc on x in this concentration range. 

To get an idea of the behavior of the average kinematic exchange A and, hence of the 
Curie temperature as a function of the Mn content and the position of the Fermi energy, 
i.e., electron concentration in the impurity band, we can make simple estimates. Consider 
the case when sp < Ecfr, then K = 1 and according to Eq. ()A23jl the levels of a pair of 
impurities separated by the distance R are split by 

fe, = ±m) = ±1^^ (29) 

where the off-diagonal lattice sum can be estimated as 

Here k is the localization parameter. (Although we may use the estimate n ~ h\^2mEcFR, 
we should not forget that it may be rather crude for really deep levels.) Then 

AiR) ^ Ao^e-''^ , 

where the value of Aq depends on the parameters of the system but is generally of the 
order of leV. 

The Poisson pair distribution truncated at small distances R < d^ is 

3R' jdl-R^\ „ ^ , 
=— exp <^ — , > , for i? > do 



giR) ={ R^ ' [ R^ }' . (31) 

0, for R< do 

This distribution accounts for the fact that the Mn impurities cannot come closer than 
at the distance do from each other. For a deep enough level we may assume k, ~ I /do 
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and calculate the dependence of the average kinematic exchange A on the Mn content 
and position of the Fermi energy with respect to the isolated impurity level Ecfr (see 
Fig. ini). The two latter quantities are measured in units Aq meaning that the kinematic 
exchange and, hence, the Curie temperature is maximal for any given value of x when 
the impurity band is half filled, i.e., Sp = E^pj^. At a; = 0.05 the kinematic exchange A 
may be about 0.14eV (if we assume that Aq = leV). The angular moment J = 2 for the 
(i^ state of the Mn impurity in GaN, hence Tc = A/Zc^, and we find that the value of the 
Curie temperature varies from a few hundred K for x = 0.01 and may exceed lOOOK at 
X = 0.05. 

It is worth mentioning that by choosing the localization radius so small, k"^ ~ do, we 
have actually found a lower bound for A. At larger values of the localization radius, we 
obtain a similar dependence of the exchange energy A on the Mn content and the Fermi 
energy position, however the absolute value of A may become essentially larger, e.g. ~0.3 
for ^ 2(io in the maximum (see Fig. [7j) meaning that the Curie temperature may 
become as high as several thousands K. 

We should emphasize here again that the averaging procedure used to obtain the results 
in Figs. El and El works well only when the Fermi energy is close to the impurity level E^p^ 
(zero on the Fermi energy axis). Far from this region one should apply the percolation 
approach proposed in Ref. |^. Then the effective kinematic exchange in this range 
may become smaller than shown in the figures, however the behavior does not change 
qualitatively, and certainly remains correct also quantitatively for Sp close to E^p^. 

To conclude this section, we have found that the optimum n-doping level for FM 
alignment of Mn spins corresponds approximately to half-filled impurity band, and one 
may well expect really high temperature magnetism in this case. The Curie temperature 
of Tq = 940K, reported in Ref.—, lies well within the range of values shown in Figs. |31and 
Unfortunately the lack of information on the impurity concentration and position of 
the Fermi level does not allow us to make a more detail comparison with the experiment. 



V. CONCLUDING REMARKS 

According to the general theory of magnetic interactions, the kinematic exchange be- 
tween localized impurity spins embedded in a semiconductor host with covalent chemical 
bonding^^ is a generic property of the two-impurity cluster. We constructed such a cluster 
which was based on an ab initio knowledge about the electronic structure of isolated Mn 
impurity in a III-V semiconductor. After some simplifications this approach led us to a 
microscopic Hamiltonian similar to a two-center Anderson modeL'^^JO That model was 
developed for metallic hosts, whereas we consider a semiconductor host. The model al- 
lows for analytical solutions, and the final equation for Tc with the effective coupling 
constant |TH|l should be compared with the phenomenological mean-field equationi22iS4 

Tc = + Xk 

where a is the lattice constant, Jpd is the phenomenological pd-exchange constant, Xh is 
the magnetic susceptibility of holes in the valence band. In principle, Eq. ()23|) and Eq. 
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FIG. 6: Dependence of the kinematic exchange, measured in units of Aq on the Mn content (x), 
and on the position of the Fermi energy relative to the level Ecfr of an isolated Mn impurity. 
We took = do. 

flH2|l consist of similar ingredients. It is natural to assume that Jpd in the phenomenolog- 
ical model stems from the hybridization. Then it may be derived in the Anderson-type 
model by means of the usual Schrieffer- Wolff transformation, so that Jpd ~ V^/ {ep — Ea). 
Looking at Eq. (fTH|) one can easily discern corresponding contributions, which in fact 
were proposed for dilute magnets four decades ago.-^ The second fraction in Eq. ()H2|1 is 
proportional to m*kp, i.e. to the Fermi momentum of holes in the valence band. Unlike 
the phenomenological model, which deals with localized spins and free holes near the top 
of the valence band, our model takes into account the change in the density of hole states 
(and therefore in their magnetic response) due to resonance scattering and impurity band 
formation. One of the results of this change is a more complicated dependence of the mag- 
netic coupling on the Mn content x than the linear law predicted by Eq. (j32|l . Another 
improvement of the mean-field theory taken into account in Eq. fl2H|l is a more refined de- 
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FIG. 7: The same as Fig. 6 but now for k = 2dQ. 

scription of the impurity magnetic moment. The actual moment J arises as a vector sum 
of the moments of the d-electrons and the bound p-hole, One should emphasize that at 
a high enough concentration all mean-field descriptions fail because the alloy approaches 
the instability region and the Mn distribution becomes strongly inhomogeneous. In the 
case of Gai^^Muj-As this is, apparently, the region of x > 0.07. 

Next, one should discuss the relation between our approach and the so called "ab initio" 
numerical calculations by means of the local-density functional method.-^^iffi Doping by Mn 
atoms is modelled in these calculations by means of a finite-size cluster of a GaAs or GaN 
host with one or several atoms replaced by Mn. Then periodic structures are constructed 
from magnetic clusters ("supercells") and the electronic bands in these artificial objects 
are calculated in the local-spin-density approximations (LSDA). So, the starting point in 
our approach and the LSDA calculations is the same. No wonder that the positions of 
CFR and DBH levels in these calculations are in good agreement with those used in our 
model, based on previous single impurity calculations . However, a direct comparison 
of the two procedures, as far as the magnetic properties are concerned, is somewhat 
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problematic. The self-consistent LSDA method results by construction only in a Stoner- 
like itinerant magnetism. Therefore, magnetic states are discussed in Refs. l29l]3fll in terms 
of exchange splitting and majority /minority spin subbands. It is difficult to discern the 
genuine kinematic double exchange in this type of calculations. 

To overcome this difficulty, Sanvito et al— tried to fit their LDSA calculation by a 
free-electron model with an effective hole mass m* and uniformly distributed impurities 
described by a model square potential containing spin-independent and exchange com- 
ponents (W{r) and respectively). Such separation is purely phenomenological. It 
ignores the resonance origin of the exchange potential. Besides, short-range impurity 
scattering cannot be described in the effective mass approximation.^^ Nevertheless, the 
estimate of the magnetic component \ J\ = 1.05 eV is in a good agreement with the corre- 
sponding parameter of our theory zV"^ / Ecfr = 1-034 eV for the parameters used in the 
fitting of Fig. El (the value of W = 0.027 eV is irrelevant because of the above mentioned 
effective mass approximation). One should emphasize, however, that the double-exchange 
coupling in our model is determined not by this parameter, but by the integral AEpM 
(|TH)). and this difference is in fact the benchmark for the discussion of the differences 
between the indirect pd exchange of Vonsovskii-Zener type^ and the kinematic double 
exchange.-^ 

A simplified picture of the valence band structure used in our model (single heavy-hole 
band with semielliptic density of states) allowed us to describe the dependence Tc{x,p) 
using a minimal number of fitting parameters. We have seen that a good quantitative 
agreement between theory and experiment could be achieved even with this very restricted 
set of parameters. With the help of a more realistic energy band scheme (e.g. by taking 
the light hole band into account), the restrictive symmetry of the density of states would 
be removed, and the fitting procedure would become more flexible. We intentionally 
imposed such severe restrictions on the model to demonstrate its explanatory capabilities. 
The above limitations will be lifted in future studies. 

Further progress in the theoretical description of ferromagnetism in Mn-doped semicon- 
ductors is intimately related to the progress in sample preparation and characterization. 
At present more or less detailed data on Tc{x,p) in p-tjpe (Ga,Mn)As are available, 
and we succeeded in a quantitative description of these data within our model (Section 
III). We expect that the same approach is applicable to p-tjpe (Ga,Mn)P, but the scanty 
experimental data do not allow us to check this expectation. As for (Ga,Mn)N, the ac- 
cumulation of experimental data is in progress, and the most actual task is to reveal 
distinctions between p and n type magnetic alloys. Recent experimental data^ confirm 
existence of magnetic order in n-(Ga,Mn)N, although the experimental state-of-the-art 
is far from providing trustworthy Tc{x) curves for theoretical fitting. According to our 
theory, FM order is expected both in p- and n-type samples, but the latter case de- 
mands a serious modification of the theory (Section IV) and a direct extrapolation of the 
semi-empirical formula (jH^ is questionable. 

Another challenging question is the possibility of ferromagnetism in elemental semi- 
conductors and II- VI compounds doped by Mn and, maybe, other magnetic ions and 
respective modification of the theory. One can mention the recent reports about mag- 
netism in Ge doped by MnSS and ZnO doped by Mn and Fe^. In the latter case one 
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deals with a wide-gap semiconductor, where magnetic isoelectronic impurities Mn^"*" and 
Fe^^ replace Zn ions. The theory presented in Section IV seems to be applicable in this 
case. Transition metal ions in Ge are as a rule interstitial impurities sometimes involved in 
formation of complex defects^"^^^, so the existing theoretical approaches should be mod- 
ified for this case. The most important conclusion from our study is that the electronic 
structure of an isolated magnetic impurity in the host material is the key to understanding 
its behavior in concentrated magnetic alloys. 
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APPENDIX A: ENERGY SPECTRUM OF TWO COUPLED IMPURITIES 

The system of four Dyson equations for the d-electron Green functions G%^i {i, i') has 
the following form within the Hartree approximation: 



\ pp' i=l,2 / 



(Al) 



Here 

gda{e)5^j = {e-Ed- UnT'' + i5sigYi{e - /i))"^5ij (A2) 

is the bare single site d-electron Green function for the t2(J electron centered at the ion z, 
taken in the Hartree approximation. ^ is the chemical potential. The Green function for 
the band electrons can be found from the Dyson equation (a similar system of equations 
for two Anderson impurities in metal was analyzed in Ref. |50|) 

p"p"'»j 

There are also impurity off-diagonal Green functions, which read 

Gl^{e)=gUe)V;,.pl^,{e). (A4) 
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The Green function G^p, describes the spectrum of the valence band electrons modified 
by the two-impurity short-range potential scattering. It satisfies the Dyson equation 

p" 

where 

9pA^) = {e-ep + i6sign{e - (A6) 

All the above equations give linear relations between the Green functions. We may 
rewrite equation ()A5|) using shorthand notations for the Green function matrices 

G6 = A-g, (A7) 

where stands for the diagonal matrix ()A6|) . Equations ()A1|) and ()A4|) take the form 

G = B ■ go (A8) 

where 

go = 



Gb 
gd 



with being a 2x2 matrix ()A2|) . The explicit expressions for the matrices A and B can 
be readily found from the Dyson equations ()A1|) . ()A4|) . and ()A5|) . 

For further calculations we need determinants of the matrices A and B, 

Q(e) = det A = q{e)'' - W^Luis)L2iie), (A9) 

where 

q{e) = l~WLu{e) 

and 

p 

is the lattice Green function. We assume that Lu = ^22- The second determinant is 

R ^ gf{e) det B = [g^\e) - V^M^.ie)]' - V'M^,{e)M^,{e). (All) 

Here 

Mfi = E e-'^P^P'^-^^G^p, = Lu + WQ-^[Ll^q + Li2L2iq + 21^^11^12^21]. (A12) 
pp' 

M22 is obtained from Eq. ()A12|) by exchanging indices 1 and 2, and 

= Yl e-^P-^^+^P'-^^G^p, = L12 + WQ-^ [LuLnq + ^2i^iiq + WLl^L^] ■ (A13) 
pp' 

Till now the calculations were made neglecting the orbital degeneracy of the impurity d 
states and within the Hartree approximation. We can generalize these calculations using 
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the Hubbard I approximation for the ci-electron Green functions (see, e.g., Ref. 1511 ). We 
assume also that the three-fold degeneracy of the impurity t2g is not lifted and the relevant 
physical quantities do not depend on the index /x enumerating these three states. The 
algebraic structure of the Dyson equation is still the same as in ()A1|) . Then considering 
the interaction of two three-fold degenerate states belonging to the two impurities. This 
leads to the 6x6 matrix 

/ a b b b\ 



a b b b 

a b b b 

b b b a 

b b b a 

\b b b a J 



(A14) 



where 



a = g-,\e) - V^KM^,{e), b = VKM^.ie) 

and K = n^s -|- . n^b is the probability that the impurity d-shell is in the nondegenerate 
state, whereas 71^^ is the probability that the impurity d-shell is in one of the three 
degenerate states. Calculating now the determinant of the matrix ()A14|) one gets the 
equation 



(A15) 



which should be used instead of Eq. ()A11|) . 

The occupation numbers n^s and nd4 for the Hubbard-like states obey a non-Fermi 
statistics, whose specific form in the case considered here is 



nd5 



CFR 



fl) 



ndi 



1 - fjEcFR - /i) 

3-2/(EcFK-/i) 



(A16) 



If the chemical potential lies below the impurity level Eqfr of the fifth electron in the 
d-shell then the Fermi distribution f{EcFR — fJ') is zero at low temperatures and ra^s = 
0, nd4: = 1/3, meaning that K = 1/3. If Eqfr < then Ud^ = 1, Ud^ = 0, and K = 1. 

The energy (0) can be found using the general property of the Green functions, con- 
necting their trace and determinant 



Tr G(e) 



d 

de 



In det G(£:) 



(A17) 



Then using Eqs. ()A9|) and ()A11|) . ()A17|) allows one to rewrite Eq. ^ in the following 
form 



= Im 



de d ^„ 



(£) + lnQ"(e)] = -Im 



de 
2^ 



[lnR"(£)+lnQ"(£)]. (A18) 



The functions R and Q depend on the combination e—isign{e—fi), hence the integration 
contour in ()A18|) can be deformed in such a way as to embrace the cut from the band 
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states and all the poles due to the localized levels with the energies below the chemical 
potential /x, i.e. occupied states. Then equation (jA18|l transforms into (jH)). 

Next we consider a property of Eq. ()A18|) . which will simplify the calculation of the 
energy and provides a better intuition to the results. Our model includes all the levels 
belonging to the valence band with the addition of the impurity d-levels, which interact 
with the band levels. Let us assume that at zero temperature the chemical potential lies 
higher than all these levels, meaning that they are all occupied. Then the total energy of 
the system is 

Etot = TtH. (A19) 

In order to calculate this trace, we may represent the operator H in matrix form using 
the noninteracting band states with the addition of the atomic ci-functions as the basis. 
Then the hybridization matrix elements will appear only in the off-diagonal positions of 
this matrix, which do not influence the result. The conclusion is that the energy Etot does 
not depend on the value of the hybridization. 

The potential scattering, Wpp, contributes to the diagonal elements of the Hamiltonian 
and may influence the value of Etot- However, we are interested here only in the indirect 
exchange between the impurities. It can be found if we consider the energy AE{Rij) and 
subtract from it the energy corresponding to two noninteracting impurities. 



1 f 
AEex = Im / de 



Ro(^) Qo(^) 



AEiocie < fi). (A20) 



where Qq{£) and Ro(£:) are obtained from ()A9|1 and ()A11|) under the assumption that 
Li2 = L21 = 0. AEioc{£ < jj) is the corresponding change of the energies of the occupied 
localized levels. Now we may conclude that the sum ()A20|) over all occupied states is 
equal to the same sum over all empty states, but with the opposite sign. Hence, 



1 

AEex = — Im / de 

TT 



. Ro(^) Qo(^)J 



AEioc{e > /i). (A21) 



Here AEioc{£ > /x) includes the empty localized levels (if any) lying above the chemical 
potential. These levels appear due to the combined action of both potential (W) and reso- 
nance {V) scattering mechanisms. Formally they can be found as zeros of the determinant 

R (unj. 

To find the contribution of the localized states to the magnetic energy one should 
simply calculate the level positions modified by the effective inter-impurity exchange due 
to hopping via empty states and their occupation. We consider here two limiting cases. 

First, we estimate the contribution of CFR levels, if they happen to lie within the 
forbidden energy gap. It means that we may expand the function R(£:) close to the energy 
^CFRa of isolated CFR level determined by the equation 

E'cpj, = E, + KV'Pn (E'cfr) , (A22) 

which describes the TM d-levels renormalized by their hybridization with the hh band. 
Indirect inter-impurity hopping results in splitting of two-impurity states and a shift of 
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localized states relative to the hh band. These levels lie in the discrete part of the spec- 
trum, where the imaginary part of the Green function (jHl) Tij = 0. Neglecting potential 
scattering, we obtain the equation 

R=[e~E,- KV'P[,i6)][e ~ E, - KV^P,^,{e)] - 9K^V^P[,{e)P,\{e) = 0. (A23) 



for the two-impurity poles in the energy gap. The solution of Eq. ()A23|) is looked for in 
the form Ecfr = E^pj^ + 6Ecfr- Now we expand the function R{Eqpj^ + 5Ecfr) up to 
the second order terms with respect to 5Ecfr and arrive at Eq. (fTT|) . 

Second, we consider the case when the DBH levels lie in the forbidden energy gap. 
Then the potential scattering W is the leading cause of the creation of the deep level. 
The energy of an isolated DBH level corresponds to a zero of the function 

q{Elj,H) = l-WPn{EUH) = ^- 

Then we look for zeros of the function R{e) at the energy Ecbh = -E'dbh + ^Edbh- 
Accounting for the fact that both functions q and Q are small in the vicinity of the energy 
Edbh the equation R = can be approximately rewritten in the form 

|q2 - W'P^, - -^W[P^,q + P^,q + 2WP,,PU | = (A24) 

{Pi2 + W[2P,,P,,q + WPU] 

with AE = E^Qjj — Ed — V^Pu. All the functions Pij are now calculated at e = E^^^. 
We first neglect the r.h.s term in Eq. ()A24|) and solve the quadratic equation 



/^From here we obtain Eq. (jlBjl for the energy shifts due to hopping between the two 
degenerate DBH levels. 

When both these levels are empty we obtain a contribution to the kinematic exchange 
by summing these two energies, extracting from them the part due to the hybridization 
with the impurity d-states, and changing the sign in the hole representation, 

^EnsH,e. = . (A26) 

Accounting for the rhs of Eq. ()A24|) will result in higher order corrections, which can be 
neglected. 
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